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Quantum Anomalous Hall Effect of Magnetic Topological Insulator Thin Films by 
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We generalize the topological phase diagram of magnetic topological insulator (TI) thin films 
in an extended parameter space comprising out-of-plane (OP) and in-plane (IP) exchange fields 
in the presence of structural inversion asymmetry (SIA) while taking a generic orbital-dependent 
spin coupling allowed for TIs into consideration. The results show that an IP held substantially 
deforms phase boundaries and generically induces the quantum anomalous Hall (QAH) effect. For 
symmetric spin coupling, extremely weak OP and IP exchange helds create the QAH state by tuning 
SIA with a gate bias. For antisymmetric coupling, the QAH phase is absent without a strong enough 
IP held. These hndings demonstrate that in the thin-hlm regime, engineering the phase boundary 
is a key process to efhciently realize and manipulate the QAH effect for nondissipative electronic 
applications. 


I. INTRODUCTION 

A thin-film topological insulator (TI) provides a planar 
platform that supports a variety of topological states and 
is suitable for exploiting related topological effects in ac¬ 
tual devices. In a thin-film configuration, surfaces states 
of a three-dimensional ( 3 D) lying on the opposite 

surfaces hybridize to create a tunneling gap. Then, the 
system alternates between the topologically-trivial nor¬ 
mal insulator (NI) phase and the nontrivial quantum spin 
Hall (QSH) phase, depending on the film thicknessA— 
The QSH state consists of two copies of quantum Hall 
(QH) states with opposite Chern numbers, arranged so 
that time-reversal symmetry is preservedJ^^— In a fi¬ 
nite system with open boundaries, there exists on each 
edge a Kramers pair of chiral edge states inside the gap. 
The helical edge mode is essentially robust against time- 
reversal-invariant nonmagnetic perturbations. However, 
magnetic disorder due to random magnetic impurities is 
detrimental to the QSH transport since they give rise to 
intraedge backscattering^i^ 

This vulnerability can be eliminated by incorporating 
ferromagnetism into the system by magnetic doping 
A strong enough exchange field normal to the surface 
drives the NI or QSH phase into the quantum anomalous 
Hall (QAH) phase, where the QH effect is activated in ei¬ 
ther of two subsystems and suppressed in the other 
The QAH phase hosts a single chiral edge mode on each 
edge, for which backscattering is completely forbidden. 
Hence, the QAH insulator exhibits perfect conduction 
that is insensitive to magnetic and nonmagnetic disorder 
as long as the gap remains intact The QAH effect is 
experimentally observed in the Bi2Se3 family of materi- 
aig 27 -_ 30 ^ and is theoretically proposed to occur in other 
related systems^^— . 

There are extensive studies in this context, where an 
out-of-plane (OP) and orbital-independent exchange field 
constitutes a main control parameter to induce the QAH 
effect for a given hybridization gap. Here, we generalize 
the topological phase diagram of magnetic TI thin films 
by considering an in-plane (IP) exchange field in addition 


to the OP field, structural inversion asymmetry (SIA) 
due to an OP electric field, and generic orbital-dependent 
exchange coupling for the Dirac model describing 3 DTIs. 
It is shown in the extended parameter space that an IP 
exchange field generically induces the QAH effect by ap¬ 
preciably deforming phase boundaries, leading to an un¬ 
conventional scheme to control the topological state of a 
thin-film TI. 


II. THEORY 

The 4 x 4 Dirac Hamiltonian describing nonmagnetic 
3 DTIs is given in momentum space by 

"Ho = £'(k) + A(k) • Txcr-I-H(k)r^, ( 1 ) 

where D = Dikz"^ D2P, A = {A2kx, A2ky, Aikz), 
B — Bq — Bikz^ — B2k'^, and k^ = kx^ + The pa¬ 

rameters in these definitions are estimated for the Bi2Se3 
familyjii^ and with fj, = x,y, z represent the Pauli 
matrices in spin and orbital spaces, respectively. To for¬ 
mulate the exchange interaction with magnetic moments 
in a homogeneous phase, we consider the generic form 
of Zeeman-type coupling allowed in the Dirac Hamilto- 
nian^i^S 

Hm = {m~\-nTz) ■ (T. ( 2 ) 

The exchange field (in energy units) is composed of the 
symmetric part m and the antisymmetric part n, by 
which the possible orbital dependence of spin coupling 
is fully describable. Thus, the total Hamiltonian consists 
of 'H = "Ho + "Hm- See Appendix A for more details. We 
first analyze the contribution due to the symmetric field 
m, which accords with the conventional approach used 
in the study of magnetic TIs. The effect of the antisym¬ 
metric counterpart n will be addressed later. 

Diagonalizing Ho under an open boundary condi¬ 
tion and projecting H onto the subspace of surface 
stateswe obtain the effective low-energy 
Hamiltonian in a thin-film configuration. The result is 



summarized as H = Hq + Hm with 


to be formulated as 
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Ho = e(k) + 752(0- X k)^ + t(k)sa; + m^cr^, (3) 
Hm = {I + cr, ( 4 ) 

in the basis constituted of top and bottom surface states 
Here, denotes the Pauli matrix in or¬ 
bital space spanned by {t,b}. Hm describes the ex¬ 
change interaction associated with the IP field m|| = 
{mx,my,Q). The contribution due to the OP field ruz 
is incorporated into Hq for convenience. Ho preserves 
rotational symmetry while it is broken by Hm- The pa¬ 
rameters are given by £ = i?o + Dk"^, t = A — Bk^, 
D = D 2 - IB 2 , B = JB 2 , I = ^ and 

J = Re (t| Tz |6). By definition, |J| reflects intersurface 
mixing, and is essentially augmented with a reduction in 
the him thickness Lz- For instance, (/, J) = (0.062,-0.7) 
at Lz = 2.3nm and (0.12,-0.35) at 3nm for Bi 2 Se 3 p^ 
whereas I — Di/Bi = 0.13 and J = 0 in the thick-hlm 
limit Lz —>■ This suggests that in the thin-him 

regime, the J term is dominant in Hm- 

The orthogonal transformation \ipi, 2 ) = (|i) ± 
arranges Ho into a block diagonal form. The eigenen- 
ergies of two blocks labeled by p = ±1 are written as 

=e + rE^P^ and E^p^ = -f'^k'^ + {t + pm-zf, where 

r = ±1 refer to the conduction and valence bands, re¬ 
spectively. The band gap Eg = 2 min | A ± uiz \ closes 
at the P point (fc = 0) when ruz^ = A^. On the other 
hand, it is shown from the Chern numbers of the asso¬ 
ciated four bands that block p becomes nontrivial when 
(A -I- pmz)B > 0,Ii^ d^d^i^3,24 ^ result, the system 

enters into the QAH phase for ruz^ > A^. If this crite¬ 
rion is not fulhlled, the system lies in the QSH phase for 
AB > 0 or the NI phase for AB < 0. These three phases 
are topologically distinct in terms of the Chern numbers 
and cannot smoothly connect without closing the gap. 

The total Hamiltonian H = Ho + Hm can be also ana¬ 
lytically diagonalized for / = 0. Defining the orthogonal 
two wavevectors ke = m||(m|| • k)/TO||^ and kt = k — k^, 
the gap function is expressed as 

E^P'> = ^J{T + pMf + 72(fc2 - jyfc/), (5) 
where T = \/+ rjj'^ke^, M = ^mz^ + and 

J^m ^ 

1 = m found from Eq. ([5]) that the gap 

closes when T — M = {1 — r])k( = kt = 0- These coupled 
conditions reduce to fc = 0 and = A^ for ruz 0, and 
predict that the QAH phase emerges in the region where 
> A^. The system otherwise belongs to the QSH 
phase for AB > 0 or the NI phase for AB < 0. Thus, 
adding Hm to Ho is formally equivalent to replacing rriz^ 
by in terms of gap closing, demonstrating that m-z 
and m\\ equally drive the phase transition into the QAH 
state. If TOz = 0, the gap remains closed as long as 
> A^, and then the system becomes metallic. 
From these considerations, we expect the Chern number 


C = - sgn(TOz)6»(mz^ -f - A^). (6) 

It is worth noting that two QAH phases classified by 
C = ±1 touch for > A^ at ruz = 0 where the 

gap vanishes. In the QSH phase, the bulk gap is open 
even for TO|| ^0, although helical edge states formed in 
a finite system are gapped by intraedge mixing. In this 
sense, the system under a small IP field is no longer in 
the ideal QSH state. In the following, we nevertheless 
refer to it as a QSH system to avoid unnecessary and 
probably confusing nomenclature. For more details, see 
Appendix B. 

A similar phase transition arising from an IP field has 
been investigated previously for paramagnetic HgMnTe 
quantum wells although its physical mechanism re¬ 
mains unresolved. To understand the underlying physics, 
we here analyze the effective Hamiltonian for a weak 
IP field. The relevant Hamiltonian is formulated as 

H = Ho + Hm{Ei^'^ — Ho)~^Hm for an unperturbed 
iv) 

eigenstate of energy L) following the second-order per¬ 
turbation formalism. H is block diagonal in the basis 
chosen for Ho- An identity term in each block is negli¬ 
gible since it has no effect on the band topology. The 
leading terms in two blocks construct the rotationally- 
invariant effective Hamiltonian, expressed as 

if = 7Sz(<T X k)z-l-t(k)sx+mzCrz, (7) 

yielding the renormalized parameters to second order 
in fc; i = A — Bk"^, A = A — J^TO||^/4A, B = 
B{1 J^m||^/4A^), and rhz = mz + /Amz- It 

is easily shown from these formulae that the original 
QAH criterion ruz^ > A^ substantially changes into 
mz^-l-J^m||^ > A^ to second order in TO| I. The renormal¬ 
ization thus reproduces the conclusion drawn from gap 
closing. 

It is also instructive to consider a generalized 3D 
parameter space (A,mz,TO||). Scanning a 2D section 
(A,TOz) along TO||, the renormalization due to TO|j is 
viewed as a continuous deformation of the QAH phase 
boundary, which consists of an equilateral hyperbola and 
a line segment joining two vertices at (±Jto||,0). As 
TO|| increases, the QAH effect is activated in the other¬ 
wise inactive region by deforming the phase boundary 
outward, i.e., TO|| intrinsically extends the QAH phase. 
However, the renormalization is no longer valid in the 
thick-film limit where A = J = 0. In this limit, top 
and bottom surface states are completely decoupled, and 
m|| merely shifts each Dirac dispersion in the opposite 
direction by an amount S 2 /(m|| x z)/ 7 . This is a clear 
distinction showing that the renormalization and the con¬ 
sequent phase-boundary deformation are characteristic of 
thin-film TIs. 
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FIG. 1. (Color online) [(a) and (b)] Band gap Eg and [(c) and 
(d)] Chern number C numerically computed as a function of 
symmetric exchange fields and In the calculation, 

SIA is neglected (V^ = 0). Solid lines represent the phase 
boundaries by which the QAH phase is separated from the 
NI phase for [(a) and (c)] = 2.3nm and from the QSH 

phase for [(b) and (d)] Lz = 3nm. 


III. NUMERICAL CALCULATION 

To examine quantitative aspects, we have performed 
numerical calculation of the Chern number using the 
Kubo formula C = -^ Shgocc /bz for a periodic 

system;^ where the summation runs over all occupied 
states at half filling, and the integration is over the Bril- 
louin zone. The momentum-space Berry curvature is 
given by 


^ ^n\^) 


m^n 


where H In) = En\n). For simplicity, we set m|| = 
(TOa;,0, 0) without loss of generality because Hq is ro- 
tationally invariant. The film thicknesses were chosen as 
Lz = 2.3nm and 3nm, for which equally |A| = 20.2meV. 
In the absence of exchange coupling, the 2.3nm- and 3nm- 
thick systems are in the NI phase {AB < 0) and in the 
QSH phase {AB > 0), respectively. Further details are 
described in our previous paper In the calculation, the 
numerical values of I and J evaluated from surface wave- 
functions were adopted. 

The numerical results are shown in Fig. |T] As ex¬ 
pected, the gap-closing line along ruz^ + 
for < A^ and ruz = 0 for > A^ separates 

the QAH phase featured by C = — sgii{mz) from the 
NI or QSH phases carrying (7 = 0. In the region where 


TTiz'^ < A^, a strong enough m|| brings about the topo¬ 
logical phase transition into the QAH phase. The for¬ 
mation of the QAH state is independently verified from 
chiral edge states in an open-boundary system and the re¬ 
sulting quantization of Hall resistance (not shown). The 
system stays in the QAH phase for > A^ irrespec¬ 
tive of rUx- In this regime, the band gap continuously 
increases with increasing nix- The QAH state is topo¬ 
logically protected and essentially immune to external 
perturbations unless the gap collapses. In this sense, the 
QAH effect is stabilized by adding an IP field through 
the gap augmentation. 

Next, we will take SIA stemming from an OP elec¬ 
tric field into consideration. The relevant Hamiltonian is 
simply given 


Hx = Vsz. ( 9 ) 

In the absence of Hm, the total Hamiltonian H = 
Hq + Hy is easily diagonalized. The gap function is 

given by [T + pM)'^ -|- (1 — where T = 

+ U^, M = y/niz"^ + and rj = It is 

found from the gap-closing condition T— M = {l — rj)k = 
0 that the QAH phase emerges when > A^ -|- U^, 
indicating that SIA tends to suppress the QAH effect 
This can be explicitly shown in (A,TOz) space from the 
phase-boundary deformation due to V. The QSH phase 
exists for AB > 0 within the interval < 

A 2 _|_ Here, Vp = l 7 l^ 0 j and fcp = yjAjB is the 
critical momentum at which t] = The NI phase 

appears otherwise, i.e., — Vp^ for AB > 0 and 

ruz^ < A^ + for AB < 0. 

Figure [2] summarizes the numerical results for Lz = 
2.3nm in the {V,mz) plane for various rux values. Solid 
lines in the figure indicate the points where the gap van¬ 
ishes. As is clearly seen, the QAH phase gradually ex¬ 
tends with increasing rUx-, and eventually covers the en¬ 
tire parameter space. A similar behavior is also observed 
for Lz = 3nm in Fig. [31 It is confirmed from these ob¬ 
servations that even in the presence of SIA, an IP field 
largely deforms phase boundaries and induces the QAH 
effect independently of the initial topology. 

A remarkable feature captured in Fig. [3] is penetra¬ 
tion of the QAH phase through the original QSH-NI 
boundary in the presence of a small m^,. To elucidate 
this peculiar behavior more quantitatively, the phase di¬ 
agram in {V,mx) space at niz = lOmeV (< |A|) is dis¬ 
played in Fig. m In this diagram, the phase bound¬ 
ary mx{V) amounts to the critical field to generate the 
QAH effect. For Lz = 2.3nm, mx{V) monotonically in¬ 
creases as a function of V from the finite lower bound 
given by mx{0) = y/ (A^ — mz^)l On the other 
hand, mx{V) decreases along the QSH-QAH boundary 
and subsequently increases along the NI-QAH boundary 
for Lz = 3nm. Importantly, mx{V) vanishes at the in¬ 
termediate critical point where these three phases con¬ 
tact. In terms of the analytical results described above, 
such a triple point appears for AB > 0 at m|| =0 and 
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FIG. 2. (Color online) Numerical values of Chern number for 
L 2 = 2.3nm as a function of V and at (a) = 0, (b) 

20meV, (c) 40meV, and (d) SOmeV. Solid lines represent the 
phase boundaries where the gap closes. 


FIG. 3. (Color online) Numerical values of Chern number for 
Lz = 3nm as a function of V and at (a) = 0, (b) 

20meV, (c) 40meV, and (d) 60meV. Solid lines represent the 
phase boundaries where the gap closes. 


= Vo^ + indicating that the QAH effect is in¬ 
ducible by extremely small and TO|| in a thin-film 
geometry. This property presents a striking contrast to 
the prerequisite niz^ > + V'^ for = which 

implies that a larger gap |A| due to a stronger intersur¬ 
face mixing obstructs the phase transition into the QAH 
state. The difficulty incidental to thickness reduction is 
circumvented by introducing a weak IP field. It is also 
noticeable that around the triple point, the topological 
state of a thin-film TI is very sensitive to V and is there¬ 
fore tunable by a moderate gate bias. The tunablity is 
intrinsic to the thin-film regime and is not easy to realize 
in thick films. 


IV. ANTISYMMETRIC EXCHANGE 
COUPLING 

Finally, we examine the effects of antisymmetric ex¬ 
change interaction. Spin coupling of this type adds the 
term 

Hn = {I + Jsj;)n:,a:,+nii- a, (10) 

to the thin-film Hamiltonian. The antisymmetric fields 
Inz and n|| can be absorbed into definitions of the sym¬ 
metric counterparts and .fin||, respectively. The re¬ 
maining term related to Jn^ is a new contribution. Sim¬ 
ilar to ruz + luz, Juz breaks time-reversal and reflection 
symmetries. The symmetry breaking is necessary for the 
QAH effecti^ Nevertheless, Juz primarily impedes the 
QAH effect. This can be explicitly shown by assuming 


c c 

- 1.0 0.0 - 1.0 0.0 



V (meV) V (meV) 


FIG. 4. (Golor online) Numerical values of Chern number at 
niz = lOmeV as a function of V and rrix for (a) = 2.3nm 

and (b) Lz = 3nm. Solid lines represent the phase boundaries 
where the gap closes. 


m|| = n|| = 0. In this simplified situation, Juz gives rise 
to only an energy splitting between two blocks, = 

Eq +pJnz, and consequently the system becomes metal¬ 
lic when Juz is sufficiently largej^ The system remains 
in an insulating state for | Jrizj < max (|A|, \mz + dnz|), 
whereas band inversion (gap closing and reopening) oc¬ 
curs in either of the two blocks when |A| < \mz + Inz\. 
Combining these two conditions, we reach the effective 
QAH criterion iJn^j < \mz + Inz\, which conflicts with 
the relation |J| > J intrinsic to thin films if rUz = 0. 
In thick films, I > | J| ~ 0 so that the surface QAH ef¬ 
fect is realizable even for niz = 0. However, edge states 
are no longer robust to external perturbations because of 
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c c 

- 1.0 0.0 - 1.0 0.0 



netic TI materials ) an IP exchange field can be 
generated by applying an external magnetic field i3|| 
parallel to the film, which rotates the magnetization 
by an angle 9 = where Bk denotes the 

anisotropy field. The auxiliary external field is suit¬ 
able to separate exchange or Zeeman coupling from or¬ 
bital effects in experiments. Moreover, in terms of the 
anisotropy energy to second order E = Ki sin^ 9 + 
K 2 sin"^ 9, the magnetization canted by an angle 9 = 
sin“^ —K 1 I 2 K 2 to the easy axis may form a ground 
state and stabilize with no external field by virtue of 
shape or surface anisotropy in a thin-film geometryi ^^do 


FIG. 5. (Color online) Numerical values of Chern number C 
for Lz = 3nm in the case of antisymmetric exchange coupling: 
(a) C{nx,nz) at P = 0 and (b) C{V,nx) at riz = lOmeV. 
Solid lines represent the phase boundaries where the gap 
closes. 

mixing with nonchiral side-surface modes;^ 

The disappearance of the QAH effect in the thin-film 
regime is reversed by introducing an IP field n\\, as shown 
in Fig.[S] Fundamentally, this effect is accountable for by 
the block mixing due to n||, which substantially changes 
the QAH criterion into J'^Uz^ + n||^ > for m = 0. 
As demonstrated above, SIA tends to quench the QAH 
effect for symmetric coupling. This is also valid for an¬ 
tisymmetric coupling. In this case, it is shown from the 
analytical diagonalization ior V 0 and =0 that 
the QAH phase is absent; the system lies in either the 
QSH phase for Vb^ > -I- J'^riz^ and AB > 0 or the NI 

phase otherwise. However, a sufficiently large n|| regener¬ 
ates the QAH effect in the presence of SIA, although the 
triple point is not located around n\\ = 0 for a small Uz- 
Thus, a nonzero IP field is essential to realize the QAH 
effect, particularly for antisymmetric coupling. This ob¬ 
servation is important in realistic situations since \nz\ is 
estimated to be comparable to, or larger than, \mz\ in a 
certain class of magnetic TIsi^ 


Appendix A: Bulk Hamiltonian 

Here, we compare in more detail the bulk Hamiltoni¬ 
ans, 'Ho and Hm, to those presented in Refs. I and 2. 
Ho is identical to Eq. (I) in Ref. I, as well as to Eq. 
(19) in Ref. 2 to second order in k. They are equally 
represented as 

Ho = H(k) -|- AiY^kz + A2{Tikx + Y2ky) ~\- H(k)F5, 

where the F matrices are defined as Fi = F 2 = 

TxCTy, F 3 = Txaz, r 4 = Ty and F 5 = r^. Ho is unitary 
equivalent to Eq. (16) in Ref. 2, expressed as 

H'o = D{k) + AiT^kz + A 2 {riky - r 2 kx) + B{k)T5, 

as explicitly shown by H'o = Hi^HoHi with Ui defined 
by Eq. (18) in Ref. 2. 

On the other hand, applying this transformation to 
Hm = (m -I- riTz) • <T yields 

H'm = {mz + nzTz)(Tz + (n|| -|- nijjTz) • cr. 

Thus, a symmetric (antisymmetric) IP field in this paper 
corresponds to an antisymmetric (symmetric) one for the 
basis chosen in Ref. 2, while symmetric and antisymmet¬ 
ric OP fields are not switched by this transformation. 


V. SUMMARY 

In summary, we have generalized the topological phase 
diagram of magnetic TI thin films by considering IP and 
OP exchange fields, SIA induced by an OP electric field, 
and symmetric and antisymmetric spin couplings with 
magnetization, on the basis of the Chern number com¬ 
puted in a multidimensional parameter space. The re¬ 
sults demonstrate that the phase-boundary deformation 
due to an IP exchange field generically induces the QAH 
effect. For symmetric coupling, the QAH state is created 
by minimal OP and IP exchange fields in the proximity of 
the QSH-QAH-NI triple point. For antisymmetric cou¬ 
pling, the QAH phase is absent without a strong enough 
IP field. 

Even when magnetocrystalline anisotropy is perpen¬ 
dicular to the film, as observed for a large class of mag- 


Appendix B: QSH phase 

In the absence of TO||, the QSH state in thin-film 
TIs is simply defined by the opposite Chern numbers 
of two blocks classified by p = ± 1 , which correspond 
to the eigenvalues of H = Since [P,Ho] = 0, p 

is a conserved quantity!^ However, it should be noted 
that the QSH state is not protected by the H-symmetry 
PHqP~^ = Hq. This is obvious from the phase tran¬ 
sition into the QAH state due to a strong enough OP 
field. 

As clearly shown in Fig. 1 (b), the bulk gap remains 
open in the QSH phase even for m\\ ^ 0. The observation 
implies the nontrivial bulk topology persists, since the 
Chern numbers cannot change for two insulating states 
adiabatically connected without closing the gap. Ana¬ 
lytically, this can be shown from the QSH criteria in the 
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c, c, 

- 1.0 0.0 10 - 1.0 0.0 1.0 


(meV) (meV) 

FIG. 6. (Color online) Nnmerical values of Chern numbers 
(a) Cl and (b) C 2 for two valence bands < £’ 2 ). Solid 
lines represent the phase boundaries where the gap closes. 
Parameters are identical to those of Fig. 1 (d), where C = 
Cl + C 2 is shown. 


FIG. 7. (Color online) Energy dispersions computed in a strip 
geometry for (a) = 0, (b) rrix = m = 20meV, (c) mj;{y) = 

m6{^ — w — \y\) and (d) mx{y) = m6{^—w — y). Parameters 
are Lz = 3nm, Ly = 2000nm, ruy = rriz = 0 and w = 30nm. 


presence of a small m||. Following the effective Hamil¬ 


tonian, Eq. (7), the system belongs to the QSH phase 
as long as ml < and AB > 0. These criteria are 
entirely verified in the QSH phase shown in Fig. 1 (d). 
Based on Eqs. (3), (4) and (8), the opposite Chern num¬ 
bers Cl = —C 2 are numerically confirmed for two valence 
bands {Ei < E 2 ) in the QSH phase, as shown in Fig. 6. 
Note that gap closing and reopening change the Chern 
number of the associated higher band. 

However, the resulting QSH effect is rather compli¬ 
cated in a finite system with open boundaries. Figure 
7 shows the energy spectra computed in a strip geome¬ 
try of width Ly with a periodic boundary condition in 
the X direction and an open boundary condition in the 
y direction. In the calculation, parameters are chosen 
to be Lz = 3nm, Ly = 2000nm and my = mz = 0. 
In the finite system, gapless helical edge modes exist for 
mx = 0 [Fig. 7 (a)], while the edge modes are gapped by 
mx = m = 20meV [Fig. 7 (b)]. 

To examine the gap-opening mechanism, we have per¬ 
formed the additional calculations that assume a nonuni¬ 
form field mxiy) partially vanishing in narrow regions 
adjacent to the boundaries at y = ±Lyj2. The width of 
the held-free region is set at w = 30nm, which is slightly 
larger than the exponential decay length of edge-state 
wavefunction £ = 26nm for mx = Note that interedge 
interaction is negligible since £ << Ly. Since w << Ly, 
the following two results are naturally expected, (i) The 
edge gap is sustained even for mx{y) if it reflects the bulk 
topology, (ii) The edge gap is suppressed for mx{y) if it 
arises from local intraedge interaction. The numerical re¬ 
sults support the latter. The edge gap almost disappears 
for mx{y) = m0{^ — w— \y\) [Fig. 7 (c)]. Moreover, for 
"mxiy) = m9{^ —w — y), the lower edge mode subjected 
to mx 7 ^: 0 is gapped while the upper edge mode little 
affected by mx tends to be gapless [Fig, 7 (d)]. These 
results confirm that in the QSH phase, intraedge mix¬ 
ing due to m|| creates the edge gap in a finite system, 
whereas the nontrivial bulk topology survives. 
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